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The SuperKamiokande group assert that they have found an oscillatory signature in atmospheric 
neutrinos through the analysis of Fully Contained Events and Partially Contained Events. We have 
performed an L/E (length/energy) analysis of Upward Through-Going Muon Events and Stopping 
Muon Events in a numerical computer simulations both with and without neutrino oscillations but 
were unable to find an oscillatory signature. We give likely explanations for the absence of the 
oscillatory signature in our simulations and its apparent presence in the SuperKamiokande data. 
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The SuperKamiokande group (hereafter SK) have re- 
ported the presence of an oscillatory signature in atmo- 
spheric neutrino data through an L/E^ analysis of neu- 
trino events occurring inside the detector, using Fully 
Contained Events and Partially Contained Events with 
energies ranging from several hundred MeV to several 
GeV [1]. Recently, it has become recognized that the 
scattering angle of the emitted lepton greatly influences 
the estimation of the direction of the incident neutrino, 
which is directly connected with the determination of L 
in the energy region of several hundred MeV to several 
GeV [2]. In the SK analysis, it is necessary for L to be de- 
cided more accurately, because the oscillatory signature 
is strongly sensitive to L/E v . 

For the Lj ' E v analysis in the SK experiment, it is more 
appropriate to analyze neutrino events occurring out- 
side the detector, such as Upward Through-Going Muon 
Events and Stopping Muon Events whose energies range 
from several GeV to several hundreds of GeV, because 
the scattering angle of the emitted lepton in the neutrino 
reactions can be neglected at such higher energies and 
consequently the direction of the emitted muon can be 
approximated as that of the neutrino events, which re- 
sults in higher accuracy in the determination of L. In this 
paper, we report the result of an analysis for an oscilla- 
tory signature in Upward Through- Going Muon Events 
and Upward Stopping Muon Events in a virtual Su- 
perKamiokande experiment by numerical computer sim- 
ulations. The particles which produce such events are 
regarded exclusively as muons (or muon neutrinos) due 
to their long flight length compared to the cascade shower 
initiated by electrons (or electron neutrinos). 

In both the present analysis of the neutrino events and 
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that of SK, a careful examination on the validity of the 
technique for the Monte Carlo Method utilized is vital, 
and we therefore start by explaining the procedures of 
the Monte Carlo Method we utilize in the present report. 
Our simulation can be regarded as a Time Sequential 
Simulation, while the SK groups is a Detector Simulation. 



We start our simulation, which is schematically shown 
in Figure 1, with the atmospheric neutrino energy spec- 
trum at the opposite side of the Earth to the detector. 
We utilize Honda's spectrum [3] as the incident atmo- 
spheric neutrino energy spectrum. We adopt a maxi- 
mum energy of 1000 GeV, and therefore the maximum 
energy of the muon emitted from the neutrino inter- 
action, here, is 1000 GeV. We calculate the range fluc- 
tuation of 1000 GeV muons by the exact Monte Carlo 
method, taking into the physical processes concerned — 
bremsstrahlung, direct pair production, nuclear interac- 
tion and ionization loss — and show the result in Fig- 
ure 2. 



We conclude from Figure 2 that it is sufficient to 
consider the neutrino events for Upward Through- Going 
Muon Events and Stopping Muon Events, which are gen- 
erated in the region within 400000 g/cm 2 of the SK de- 
tector (less than 2000 meters), because neutrino inter- 
actions further than 2000 m from the detector could not 
contribute physical events into the detector. 



Then, we define Ni n t(E v ,t,cos6 v )dt, the interaction 
neutrino energy spectrum at depth t from the detector 
underground for the incident neutrino with energy E v 
from the zenith angle 0„, in the following. 
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FIG. 1: Schematic illustration of the experiment. 
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FIG. 2: Range Energy Fluctuation of 1000 GeV muons. 
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FIG. 3: Interaction energy spectrum without neutrino oscil- 
lation. 
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FIG. 4: Interaction energy spectrum with neutrino oscilla- 
tion. 
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(1) 



Ai(E u ,ti, pi) Xi(E u ,ti, pi) Xi(Ep,ti, pi) 

where \i(E v ,ti, pi) denotes the mean free path of the 
neutrino with energy E v at the distance U from the 
opposite surface of the Earth whose density is pi and 
\i(Ev,ti, pi) denotes the corresponding mean free path 
of the anti- neutrino whose energy is given by Ep, These 
mean free paths are calculated from the deep inelastic 
scattering cross sections [4]. The density of the Earth is 
adopted from the Preliminary Earth Model for the den- 
sity profile [5]. 



The distribution functions for L/E v , obtained using 
Eq.(l), are given in Eq.(2) and Eq.(3) in the cases with- 
out neutrino oscillation and with neutrino oscillation, re- 
spectively. 



N null-osc | _^ cosd \ d (L 

E v I \E V 



N v {E v ,L,cos9 v )d 



L 



(2) 



N V {E V , L, cos6„)P(^ — » v»)d ( — ) (3) 



where P{v^ — > Vu) denotes the survival probability 
for a neutrino in the presence of neutrino oscillations, 
which is given as 
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FIG. 5: Distribution of L/E^ for Upward Through-Going FIG. 6: Distribution of L/E^ for Upward Stopping Muon 
Muon Events. Events. 



= 1 - sin 2 29s 



l.27Am 2 (eV 2 )L(km) 
E„{GeV) 



(4) 



Here, we adopt sin 2 28 — 1.00 and Am 2 = 2xlO~ 3 eV 2 , 
as obtained from the SK experiment. 

In Figures 3 and 4, we give the interaction neutrino 
energy spectrum without, and with, neutrino oscillations 
as defined by Eq.(2) and Eq.(3), respectively. It is clear 
from Figure 4 that the effect of the neutrino oscillation 
for the SK parameters does not appear in the horizontal 
direction, cos9 v = 0.0, due to short path length for the 
traversed neutrino while the effect clearly appears the 
vertical case, cos9 v = 1.0 

For Upward Through- Going Muon Events and Stop- 
ping Muon Events, we can assume that the direction of 
the incident neutrino is the same as that of the emitted 
muon, because the scattering angle of the emitted muon 
can be neglected due to its high energy. The simulation 
procedures for the events concerned for a given zenith 
angle of the incident neutrinos are as follows. 

Procedure A: Using £, a uniform random number be- 
tween and 1, we sample the energy of the inci- 
dent (anti-) neutrino, which is obtained from the following 
equation: 



•/"Wmin Nint(Eu, t, cos9 v )dt 
J E "' m "° N int (E v ,t, cos9 v )dt 



(5) 



Procedure B: For £„(p), the energy of the (anti- 
neutrino obtained by Procedure A, we define 



o{E v ) 



a{E v ) + a(E p ) 



where a(E„) and a(E p ) denote the total cross-section 
of the neutrino and the anti-neutrino, respectively. Then 



we sample the random number £ again and if £ < £x we 
take the incident lepton to be a neutrino, otherwise we 
take it to be an anti-neutrino. 

Procedure C: We decide the interaction point of the 
(anti-) neutrino event determined by Procedure B in the 
range (0, 2000) meters. This is the distance from the 
detector to the interaction point and is obtained simply 
by sampling a uniform random number between (0,1) as 
this range is many orders of magnitude smaller than the 
mean free path of the neutrino concerned. 

Procedure D: We sample , the energy of the (anti- 
)muon emitted for the deep inelastic scattering by using 
the uniform random number between (0,1), which is log- 
ically same as in Eq.(5). 

Procedure E: For the (anti-)muon whose energy and 
production point is determined from Procedures C and 
D, we examine the behavior of the trajectory of the 
(anti-)muon toward the SK detector in a stochastic man- 
ner. Namely, each individual muon is pursued by taking 
into consideration bremsstrahlung, direct pair produc- 
tion, nuclear interaction, and ionization loss without uti- 
lizing the average behavior of the muon concerned. As 
the result, we determine which category each individual 
muons falls into: [a] stopping before it reaches the detec- 
tor, [b] stopping inside detector, or [c] passing through 
the detector. 

We repeat Procedures A to E and obtain the neutrino 
events concerned for a given live time for the real experi- 
ment. In our computer numerical simulation, we have ac- 
cumulated the events concerned which correspond to the 
real live time for SK. For each neutrino event, we know 
E v , the energy of the parent neutrino, E^, the energy of 
the daughter (anti-)muon, cos9 v (cos# M ), the direction of 
both the incident (anti-)neutrino and the emitted (anti- 
)muon, and L, the distance between the interaction point 
of the neutrino events and the opposite side of the Earth. 

As only E^ can be measured in the actual SK experi- 
ment, and not E Vl we give the frequency of the number 
of the events as a function of L/ ' E^ for Upward Through- 
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FIG. 7: E u vs scatter plot for Upward Through-Going 
Muon Events. 



Neutrinos and induced muons (Upward stopping muon) 
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FIG. 8: E v vs _B M scatter plot for Upward Stopping Muon 
Events. 



Going Muon in Figure 5 and the corresponding quan- 
tity for Stopping Muon Events in Figure 6 in both the 
cases without and with neutrino oscillations. It is clear 
from those figures that no oscillatory signature is appar- 
ent, with almost no difference between the cases with 
and without oscillation when plotted against the mea- 
sured SK parameter. Of course, according to the logic 
adopted by SK, the oscillatory signature should appear 
in the function of L/E v , not L/E^. 

The energy of the emitted lepton for a given energy of 
the neutrino is generated in the deep elastic scattering, 
producing a wider energy distribution. In Figures 7 and 
8, we give the scatter plots between E v and for Up- 
ward Through- Going Muon Events and Stopping Muon 
Events, respectively. It is clear from the figures that 
the fluctuations in the energy distributions are not small. 
The density of points in the figures is proportional to the 
numbers of events. 

In Figures 9 and 10, we give the corresponding distri- 
butions for EjE v to Figures 5 and 6 for L/E^. From the 
comparison of Figures 5 and 6 with Figures 9 and 10, 



Distribution of L(km)/E_neutrino(GeV), (upward through-going muon) 
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FIG. 9: Distribution of LjE v for Upward Through-Going 
Muon Events. 



Distribution of L(km)/E_neutrino(GeV). (upward stopping muon) 
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FIG. 10: Distribution of L/E v for Upward Stopping Muon 
Events. 



it can be concluded that oscillatory signatures are not 
observed in both EjE v and L/E^ distributions. Also, 
there are only small differences between the L/E distri- 
butions with and without neutrino distribution, as seen 
in Figures 5 and 6 and in Figures 9 and 10, which is 
consistent with our previous conclusion that there is no 
evidence for neutrino oscillation through the analysis of 
the zenith angle distribution for Upward Through- Going 
Muon Events and Stopping Muon Events [6]. 

In Figures 11 and 12, we give the ratios of "with 
oscillation" to "without oscillation" as the function of 
EjE v for Upward Through-Going Muon Events and Stop- 
ping Muon Events, respectively. Generally, the fact 
that fluctuations are rather larger comes from relatively 
small number of events, and further, the fluctuations are 
larger in Upward Stopping Muon Events than in Upward 
Through- Going Muon Events, as they must be. In Fig- 
ures 13 and 14, we give the corresponding quantities for 
L/E^. Comparing Figures 11 and 12 with Figures 13 
and 14, the fluctuations are larger in the latter than in 
the former, which is easily understandable when consid- 
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FIG. 11: The ratio of "oscillation" to "null oscillation" of 
L/E u distribution for Upward Through- Going Muon Events. 



FIG. 13: The ratio of "oscillation" to "null oscillation" of 
L/Ep distribution for Upward Through-Going Muon Events. 
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FIG. 12: The ratio of "oscillation" to "null oscillation" of 
L/E v distribution for Upward Stopping Muon Events. 
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FIG. 14: The ratio of "oscillation" to "null oscillation" of 
L/E^ distribution for Upward Stopping Muon Events. 



ering the energy distribution of for a definite E v . It is 
clear from these figures that it is exceedingly difficult to 
extract some positive evidence for neutrino oscillations. 
In Figures 15 and 16, we give scatter plots between L 
and E v with and without neutrino oscillation for Upward 
Through- Going Muon Events and Stopping Muon Events, 
respectively. Also, we give the corresponding plots be- 
tween L and in Figures 17 and 18. 

Here, we comment on the evidence for an oscillatory 
signature claimed by SK. In our opinion, it is practically 
impossible to observe such an oscillatory signature in at- 
mospheric neutrino oscillation for the following reasons: 

(i) The L/E distribution for the muon events oc- 
curring inside the detector in the SK experiment, 
N(L/E v )d(L/E v ) is given as, 



d(^rj J dcos9 v D nsp (E u ,cos9 v ,L(cos9 u )) 
xP(i^ — yp^) / a u ^ fJ ,(E v ,E /J ,)dE lJ 



where D nsp (E v , cos9 v , L(cos9„)) is the differential neu- 
trino incident energy spectrum at the given zenith an- 
gle, 9 V , and L, which is also the function of cos9 v . 



denotes the probability function for the 



neutrino oscillation and a^^^E^^E^jdE^ denotes the 
differential cross-section for the incident neutrino to pro- 
duce the emitted muon inside the detector. The character 
of a oscillatory signature obtained by SK comes exclu- 
sively from the part of P{Vfi — ► v^) in Eq.(6). However, 



Ply* 



in Eq.(6) is imbedded into the steep neu- 



trino energy spectra like that shown in Figure 4 and is 
itself a continuously and strongly varying function with 



(6) 
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FIG. 15: L vs E v scatter plot for Upward Through-Going 
Muon Events. 



FIG. 17: L vs E^ scatter plot for Upward Through-Going 
Muon Events. 
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FIG. 16: L vs E v scatter plot for Upward Stopping Muon 
Events. 



FIG. 18: L vs Efj, scatter plot for Upward Stopping Muon 
Events. 



L. As a result, it can not be separated from the neutrino 
energy spectrum. 

(ii) The SK Detector Simulation assumes that the di- 
rection of the incident neutrino is the same as that of 
the emitted lepton [7,8]. Thus, SK could treat neither 
the effect of the azimuthal angle of the emitted lepton 
nor the effect of the backscattering over the zenith an- 
gle of the emitted lepton and consequently could not 
determine L reliably, which is directly connected with 
N(L/E u )d(L/E v ). More concretely, under the combina- 
tion of the two fundamental parameters in the neutrino 
oscillation derived by SK, about 40% of the upward going 
muons originate from downward going neutrinos, while 
about 10% of the downward going muons originate from 
upward going neutrinos [9]. Such mutual mixing in the 
direction of the incident neutrinos may be maximized 
near L/E v — 150 ~ 500 (km/GeV), where one find the 
dip-like phenomena of the N(L/E„)d(L/E v ) in Figure 3 
in the SK paper [1]. If we determine the direction of 
the incident neutrino correctly, then the dip will disap- 



pear and N(L/E u )d(L/E u ) will show the something like 
behavior given in Figures 9 and 10. Our subsequent pa- 
per on the L/ E analysis for Fully Contained Events and 
Partially Contained Events will be published elsewhere. 
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